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To  make  the  polymer  electrolyte  membrane  fuel  cells  (PEMFC)  commercially  viable,  further  reduction 
in  cost  and  improvement  in  performance  are  required.  In  this  work,  an  innovative  design  of  a  PEMFC 
with  multiple  catalyst  layers  (CLs)  is  considered  and  the  design  variables,  weight  fraction  of  platinum  on 
carbon  ( fpt ),  platinum  loading  (mpt),  ionomer  loading  C/jon0mer)  and  thickness  of  all  the  CLs,  are  optimized 
for  cost  reduction  and  performance  enhancement.  In  the  first  optimization  study,  the  cell  performance 
is  maximized.  The  maximum  current  density  of  an  optimized  PEMFC  with  four  CLs  shows  a  significant 
improvement  over  the  base  case  design  at  all  operating  voltages.  In  another  optimization  study,  the 
platinum  loading  of  the  PEMFC  with  multiple  CLs  is  minimized.  With  an  increase  in  the  number  of  catalyst 
layers,  the  platinum  required  to  achieve  the  base  case  design  current  density  is  reduced.  Using  four  CLs,  a 
reduction  of  17%  (at  0.15  A  cm-2)  to  60%  (at  0.7  A  cnrr2)  in  platinum  loading  is  achieved  in  comparison  to 
the  base  case  design.  In  addition,  the  maximum  power  density  of  the  PEMFC  with  multiple  CLs  is  found 
to  be  superior  to  that  of  an  optimized  PEMFC  with  a  single  CL  at  all  voltages. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  cathode  catalyst  layer  (CL)  is  one  of  the  most  important  lay¬ 
ers  of  a  polymer  electrolyte  membrane  fuel  cell  (PEMFC).  Some  of 
the  issues  that  have  hindered  commercialization  of  the  PEMFCs  are: 
(i)  under-utilization  of  the  electro-catalyst  (ii)  lower  performance 
due  to  mass  transfer  losses,  ionic  and  ohmic  losses,  and  inadequate 
removal  of  liquid  water  from  the  reaction  sites.  To  mitigate  these 
problems,  we  investigated  an  innovative  design  of  a  PEMFC  with 
multiple  catalyst  layers  in  one  of  our  previous  papers  [1].  There 
are  few  other  experimental  and  numerical  studies  that  have  also 
highlighted  the  advantages  of  multiple  catalyst  layers  over  the 
conventional  single  layer  [2-4].  The  motivation  for  using  multi¬ 
ple  catalyst  layers  is  to  vary  the  design  parameters  of  a  CL  spatially 
according  to  the  operational  requirement  of  a  PEMFC.  The  CL  adja¬ 
cent  to  the  diffusion  medium  should  be  of  higher  porosity  than  the 
other  CLs  as  the  concentration  of  oxygen  decreases  in  a  CL  from  the 
diffusion  medium-CL  interface  towards  the  polymer  membrane. 
The  CL  adjacent  to  the  polymer  membrane  should  contain  more 


Abbreviations:  PEMFC,  polymer  electrolyte  membrane  fuel  cell;  GC,  gas  chan¬ 
nel;  GDL,  gas  diffusion  layer;  MPL,  microporous  layer;  CL,  catalyst  layer;  MEM, 
membrane;  MEA,  membrane  electrode  assembly;  ORR,  oxygen  reduction  reaction. 
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ionomer  than  the  other  CLs.  Furthermore,  liquid  water  should  be 
removed  without  causing  significant  mass  transport  and/or  ohmic 
losses.  Because  of  the  difficulty  in  manufacturing  such  a  continu¬ 
ously  graded  CL,  a  combination  of  layers  can  be  synthesized  where 
each  layer  is  manufactured  with  different  design  parameters.  Our 
previous  paper  showed  the  performance  of  the  PEMFC  with  multi¬ 
ple  CLs  is  superior  to  the  PEMFCs  with  a  single  CL  [1  ].  Motivated  by 
these  results  we  have  performed  optimization  studies  of  a  PEMFC 
with  multiple  CLs  to  further  improve  its  performance  and  reduce  its 
cost.  The  typical  design  parameters  of  the  CL  are:  weight  fraction  of 
platinum  on  carbon  (fpt ),  platinum  loading  ( mpt )  and  ionomer  load¬ 
ing  (Jionomer )•  For  optimizing  these  design  parameters,  a  detailed 
steady  state  model  of  a  PEMFC  cathode  with  multiple  layers  is 
required.  In  this  work,  we  have  considered  our  previously  devel¬ 
oped  model  with  few  modifications.  The  basic  PEMFC  configuration 
remains  similar  to  before.  The  model  considers  liquid  water  in 
all  the  layers.  The  catalyst  layer  microstructure  is  modeled  as  a 
network  of  spherical  agglomerates.  This  characterization  is  based 
on  a  previous  study  [5]  that  compared  macro-homogenous  and 
spherical  agglomerate  characterizations  and  observed  that  spher¬ 
ical  agglomerate  characterization  is  a  better  representation  of  the 
CL.  For  improved  water  management,  a  thin  micro-porous  layer 
is  considered  between  the  gas  diffusion  layer  (GDL)  and  the  first 
catalyst  layer. 

A  few  papers  are  available  in  the  open  literature  on  the  cathode 
CL  optimization.  Song  et  al.  [6]  have  obtained  optimal  distributions 
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Nomenclature 

ci 

concentration  of  species  i  in  region  k  (molm-3) 

Q,o 

Inlet  concentrations  of  i  (mol  m-3) 

Q)2  ns 

concentration  of  dissolved  oxygen  at  the  interface 
of  ionomer  and  agglomerate  (molm-3) 

nmem 

concentration  of  liquid  water  in  the  membrane 
(molm-3) 

Df’k 

effective  diffusivity  of  the  species  i  in  region  k 

(m2  s-1) 

Dmem 

diffusivity  of  oxygen  in  ionomer  (m2  s-1 ) 

Dr 

diffusivity  of  liquid  water  in  the  membrane  (m2  s-1 ) 

fionomer 

weight  fraction  of  ionomer  in  the  catalyst  layer 

fPt 

weight  fraction  of  platinum  on  carbon 

F 

Faraday’s  constant  (Cg-1  equiv.-1) 

1-base 

current  density  for  base  case  design  conditions 
(A  m-2  Pt-1 ) 

i-cell 

cell  current  density  (Am-2  Pt-1 ) 

iopt 

current  density  for  optimized  design  (Am-2  Pt-1 ) 

/w 

interfacial  transfer  of  water  between  liquid  and 
vapor  (molm-3  s-1) 

Jk 

local  flux  due  to  diffusion  of  species  i  in  region  k 
(molm-2  s-1) 

kc 

condensation  constant  (s-1 ) 

ky 

evaporation  constant  (atm-1  s-1 ) 

mpt 

platinum  loading  inside  the  catalyst  layer 
(kg  Ptm-2  CL-1) 

Mw 

molecular  weight  of  water  (g  mol-1 ) 

n 

number  of  electrons  taking  part  in  the  oxygen 
reduction  reaction 

Nw,k 

flux  of  liquid  water  in  region  k  (mol  m-2  s-1 ) 

Pc 

capillary  pressure  (atm) 

Pw 

partial  pressure  of  water  vapor  (atm) 

rjsat 

Vw 

saturation  pressure  of  water  vapor  (atm) 

Q 

switching  function 

ragg 

agglomerate  radius  (m) 

R 

universal  gas  constant  (J  mol-1  K-1 ) 

Ro2 

rate  of  oxygen  reduction  reaction  per  unit  volume 
of  the  catalyst  layer  (mol  m-3  s-1 ) 

sk 

liquid  water  saturation  level  in  region  k 

Tcell 

cell  temperature  (I<) 

tcL 

thickness  of  the  catalyst  layer  (m) 

Vcell 

cell  voltage  (V) 

yw 

mole  fraction  of  water  in  the  gas  phase 

Greek  letters 

Oia 

vapor  activity  in  the  gas  phase  in  the  anode  catalyst 
layer 

Oic 

vapor  activity  in  the  gas  phase  in  the  cathode  cata¬ 
lyst  layer 

£k 

void  fraction  inside  region  k 

& ionomer 

fraction  of  volume  occupied  by  the  ionomer  inside 
the  catalyst  layer 

$ 

effectiveness  factor 

Keff.c 

effective  proton  conductivity  in  the  catalyst  layer 
(mhom-1) 

Keffmem 

effective  proton  conductivity  in  the  membrane 
(mhom-1) 

iA 

e/e 

electric  conductivity  in  region  k  (S  m-1 ) 

Keffk 

Kele 

effective  electric  conductivity  in  region  k  (S  m-1 ) 

Pc 

density  of  carbon  (kgm-3) 

Pionomer 

density  of  ionomer  (kgm-3) 

PPt 

density  of  platinum  (kgm-3) 

pw  density  of  water  (kg  m-3 ) 

ifs  Thiele  modulus 

Subscripts 

i  index  for  the  species:  02,  N2,  H20 

k  index  for  the  region:  diffusion  layer,  micro-porous 

layer,  catalyst  layer 


of  platinum  and  Nation  for  the  maximization  of  current  density  at  a 
given  voltage.  The  authors  concluded  that  the  optimal  distribution 
of  Nation  content  is  a  linearly  increasing  function  and  the  optimal 
distribution  of  platinum  is  a  convex  increasing  function  across  the 
thickness  of  the  catalyst  layer.  Lin  et  al.  [7]  optimized  the  chan¬ 
nel  to  width  ratio,  porosity  of  the  GDL,  and  porosity  of  the  catalyst 
layer.  Secanell  et  al.  [8]  performed  multivariable  optimization  stud¬ 
ies  of  the  PEMFC  cathode.  The  authors  considered  platinum  loading, 
volume  fraction  of  the  ionomer  in  the  agglomerate,  fraction  of  the 
platinum  on  carbon,  and  porosity  of  the  GDL  as  decision  variables. 
The  authors  suggested  that  the  performance  at  medium  and  high 
current  densities  can  be  improved  by  increasing  the  ionomer  con¬ 
tent  and  reducing  the  carbon  and  platinum  loadings.  The  authors 
did  not  consider  the  liquid  water  effect  on  the  cell  performance. 
The  liquid  water  generated  in  the  ORR  reaction  plays  a  major  role 
especially  in  the  mass  transfer  limiting  regions.  Later  Jain  et  al.  [9] 
reformulated  the  agglomerate  model  of  Secanell  et  al.  [8]  into  a  con¬ 
densed  form  for  minimization  of  the  platinum  loading  at  various 
current  densities.  They  obtained  optimum  platinum  distribution 
along  the  CL  width  by  dividing  the  CL  into  various  zones. 

Our  group  recently  published  a  multivariable  optimization 
study  of  a  PEMFC  cathode  considering  the  effect  of  liquid  water 
[10].  In  the  current  work,  multivariable  optimization  studies  of  a 
PEMFC  with  multiple  CLs  are  performed.  Two  objective  functions 
are  considered  here.  In  the  first  study,  the  objective  is  to  maximize 
the  cell  current  density  for  a  given  voltage.  The  second  study  focuses 
on  minimization  of  the  platinum  loading  for  a  given  current  den¬ 
sity.  The  design  parameters  of  all  the  catalyst  layers  are  the  decision 
variables  for  both  the  optimization  studies. 

2.  Model  description 

Because  of  some  modifications  to  the  previously  developed 
model  [1  ]  and  to  ensure  completeness,  the  model  used  in  the  opti¬ 
mization  studies  is  briefly  described  here.  The  schematic  of  the 
model  domain  considered  in  this  study  is  shown  in  Fig.  1.  The 
cathode  flow  field  consists  of  parallel  channels.  The  reactants  fed 
to  the  channels  diffuse  to  the  CL  through  the  gas  diffusion  layer 
(GDL)  and  micro-porous  layer  (MPL).  The  two-dimensional  two- 
phase  model  with  four  CLs  considers  transport  of  the  reactant 
gases,  02,  N2,  and  H20  (v),  in  the  GDL,  MPL,  and  CLs.  In  addi¬ 
tion,  liquid  water  transport  in  the  membrane  (considering  both 
electro-osmotic  drag  and  back  diffusion),  CLs,  MPL,  and  GDL  is  con¬ 
sidered.  Furthermore,  protons  transport  in  the  membrane  and  CLs, 
and  electrons  transport  in  the  GDL,  MPL,  and  CLs  are  modeled.  Elec¬ 
trochemical  oxygen  reduction  reaction  (ORR)  is  considered  in  all 
the  CLs.  The  model  assumptions  are:  (i)  isothermal,  isobaric,  and 
steady-state  operation;  (ii)  water  generated  due  to  ORR  is  in  liq¬ 
uid  form;  (iii)  physical  properties  of  ionomer  in  all  the  CLs  and  in 
the  membrane  is  same;  (iv)  negligible  contact  resistance  between 
the  cathode  layers.  The  model  equations  and  boundary  conditions 
are  shown  in  Tables  1  and  2,  respectively.  The  system  of  equa¬ 
tions  is  solved  using  a  combination  of  MAPLE®  and  MATLAB®. 
The  partial  differential  equations  (PDEs)  along  with  the  bound¬ 
ary  conditions  are  written  in  MAPLE®.  The  PDEs  are  discretized 


M.  Srinivasarao  et  al.  /  Journal  of  Power  Sources  206  (2012)  197-203 


199 


Excess  reactants  and 
water  droplets  out 


Flow 

field 


C  C  C  C  M 
L  L  L  L  E 
1  2  3  4  M 


Reactants  gases  in 


02>  N*  H20  (v) 


Liquid  water 


O  Spherical  agglomerate 

. ^  Diffusion  due  to  cone,  gradient 

- -►  Electro  osmotic  drag 

GDL-Cas  Diffusion  Layer 
M PL-Micro-Porous  Layer 
CL  -  Catalyst  Layer 
MEM  -  Membrane 


Fig.  1.  Schematic  of  the  model  domain. 


Table  1 

Conservation  equations  (C=  concentration,^  gas  flux,  N=  liquid  flux). 


Variables 

Gas  channel 

Diffusion  layer  (DL) 

Micro  porous  layer  (MPL)  Catalyst  layers  (CLs) 

Polymer  membrane 

Co2 

-|(Co2u)-v-io2  = 

0 

-V-(— Dg'dVCo2)  =  0 

-V(-D«mVCo2) 

=  0  -V.(-D^cVCo2)- 

Kc 

>2=0  - 

Cn2 

-^(Cn2u)  -  V  ■ Jn2  = 

0 

-V.(-D«dVCN2)  =  0 

-V.(-D«mVC„2) 

=  0  -V.(-D^cVCn2)  = 

0 

- 

Ch2o 

~-§(Ch2ou)  ~  V  Jh2o 

—  iw  =  o 

-V.(-D«dVCH2o)-W 

=  0  -V-(-Df^VCH2o)-Iw  =  0  -V.(-D^VCh2o)- 

-Iw  =  0 

s 

-^r|(su)-VNw,d  +  /w  =  0  - 

-  V  ■  N^,d  +  lw  =  0 

-  V  ■  Nw,m  +Iw~0 

—V  ■  Nw,c  +  Iw  +  2  Rq2  ■ 

—  0  —  V  ■  Nw,mem  ~  0 

<\>r 

- 

- 

~ 

K«cvVr  +  nFR  02  = 

0 

_Keff,memv2(pr  =  0 

</>s 

- 

=  0 

-Kf'mv24,s  =  o 

<fV20s-nFRO2  = 

0 

Where  Iw  =  kc(ek(  1  -  s)/RTceli)yw(pw 
oxygen  consumption). 

-  Psw)q  +  kv(ekspw/Mw)(pw  -  pyx  1  - 

q )  (interfacial  transfer  of  water  between  liquid  and  vapor  phases)  and  Rq2  =  -£/<™nCo2  Ins  (rate  of 

Table  2 

Boundary  conditions  in  thickness  direction. 

Variables 

Co2 

Cn2 

Ch2o 

s 

0r 

0s 

Entrance 

LJ 

II 

CJ 

p 

II 

p 

o 

o 

II 

j? 

o 

0 

- 

- 

GC/DL 

II 

II 

o 

II 

sGC=sDL 

- 

0S“  =  Vcdl 

DL/MPL 

rDL  _  rMPL 

02  02 

y 

II 

/-DL  _  rMPL 

^h2o  h2o 

pDL^pMPL 

<f‘'V0s|DI=/cfmV0s|MPl 

og 

II 

II 

if 

II 

o 

h]DL  _  MPL 

nw  ~  JV 

- 

0sDL  =  0sMPL 

MPL/CL1 

*1 

II 

onp 

II 

if 

rMPL  _  rCL\ 

^H20  ~  ^ H20 

r,MPL  _  nCLl 
fc  —  rc 

O 

II 

> 

<f’mV0s|MPL=<faiV0siCL  1 

P 

II 

[op 

d^ 

II 

iMPL  _  rCLl 

jh2o  ~jh2o 

mMPL  _  ajCL\ 

0SMPL  =  0?1 

CL1/CL2 

U° 

II 

CJ° 

rCL\  rCL2 

^n2  u n2 

*P 
o  ^ 

II 

o 

pm  =PCL2 

0r|CLl  “  0r|CI2 

<fL1V0s,cLi=4’a2V0s|a2 

to5  P 

II 

fpp 

d^ 

II 

d^ 

Jrp 

to  ^ 

II 

[srp 
o  w 

NW=NW 

V  </)r|C£1  =  iceff<CL2  V  0r|CL2 

0sCil  =  0sCi2 

CL2/CL3 

II 

u° 

II 

/-CL2  rCL3 

^ h2o  ~  h2o 

p?2=p CL3 

0r|CL2  =  0r|CI3 

<fa2V0s|CL2=<fa3V0s|a3 

dc? 

ii 

dc? 

t?p 

II 

jrg 

Jrp 

tO[° 

II 

[3cp 
O  w 

NCL2  _  NCL3 

^Vte  =  ^V^|Ci3 

0sCI2  =  0sCL3 

CL3/CL4 

u° 

II 

u° 

II 

/— CL3  /-CL4 
^H20  —  ^ H20 

P°cL 3  =PCcL4 

0r|CL3  =  0r|CL4 

,  ej T,  CL3  w  ,  e/L  CL4  w  j. 

Kde  V0S|a3  =  Kde  V0S|CL4 

erg 

ii 

op 
10  E 

Jrp 

II 

^p 
^  E 

3°. 

05 

II 

mO 

d£" 

NCL3  =  NCL4 

^a3V0r|CL3=^CL4V^a4 

0?3  =  004 

CL4/MEM 

VC?4  =  0 
°2 

< 

^p 

II 

o 

< 

^P 
o  ^ 

II 

o 

CMEM  =  C^V) 

0r|CL4  =  0r|  MEM 

V0?4  =  0 

NCL4  _  NMEM 

Keff’CL4  V  (pr\CL4  =  V  0r|MEM 

MEM/ANODECL 

qMEM  _  ^CL,EQ^anodej 

0r  =  0 
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Table  3 

Base  case  design  decision  variables. 


Decision  variable 

Base  case  design 

mpt  (mg cm-2) 

0.4 

frt 

0.2 

fionomer 

0.34 

tCL  (|xm) 

20 

Table  4 

Range  of  design  parameters  for  optimization. 

Decision  variable  (i  =  1  —4) 

Lower  bound 

Upper  bound 

mpt,i  (mg cm-2) 

0.0125 

0.1  and  0.25 

fpt,i 

0.05 

0.95 

fionomer, i 

0.05 

0.95 

ku  (fxm) 

2.5 

5 

in  MAPLE®  using  forward  difference  technique.  Analytical  Jaco- 
bians  are  generated  and  then  the  system  of  equations  along  with 
the  analytical  Jacobians  is  exported  to  MATLAB®.  The  optimization 
is  performed  in  MATLAB®  using  the  function  “fmincon”  available 
in  the  MATLAB’s  optimization  tool  box.  The  constrained  nonlin¬ 
ear  optimization  has  been  carried  out  by  a  feasible  path  approach 
in  which  the  PEMFC  model  equations  are  converged  at  each  iter¬ 
ation  by  a  general-purpose  nonlinear  solver  “fsolve”  available  in 
MATLAB®.  “fmincon”  takes  about  20-25  iterations  to  converge 
depending  on  the  operating  voltage.  For  each  iteration,  “fsolve” 
takes  about  20  iterations  to  converge.  The  details  of  modeling  and 
optimization  framework  can  be  found  in  our  previous  publications 
[1,10]. 

3.  Optimization  studies 

The  base  case  design  parameters  are  given  in  Table  3.  The  ranges 
of  the  decision  variables  considered  in  this  study  are  provided  in 
Table  4.  The  upper  bound  on  the  thickness  of  the  ultrathin  CL  (ta) 
is  kept  at  5  |xm  so  that  the  total  thickness  of  all  the  CLs  would  not 
exceed  20  [xm.  Constraints  are  imposed  on  the  volume  fractions  of 
voids,  ionomer,  and  the  solids  in  all  the  CLs  to  ensure  that  these 
fractions  are  non-negative  and  less  than  one. 

3.1.  Problem  formulation-l 

In  formulation  I,  the  cell  current  density  is  maximized  for  a 
given  operating  voltage.  To  understand  the  effects  of  platinum 
loading  in  each  CL,  the  problem  is  solved  with  different  upper 
bounds  and  constraints.  Since  the  platinum  loading  in  each  CL  in 
the  base  case  design  is  0.1  mg  cm-2,  the  upper  bound  on  platinum 
loading  for  each  CL  is  kept  at  0.1  mg  cm-2  in  case  1.  This  ensures 
that  the  maximum  overall  platinum  loading  on  the  cathode  side  is 
0.4  mg  cm-2.  In  the  second  case,  the  upper  bound  in  each  CL  is  kept 
at  0.25  mg  cm-2.  The  overall  platinum  loading  is  1  mg  cm-2  in  this 
case.  Cases  1  and  2  will  be  referred  to  as  low  Pt  case  and  high  Pt 
case,  respectively,  in  our  subsequent  discussion. 

Objective  function:  Maximization  of  iCeii  at  a  given  VCeu 

Decision  variables  :  fptJ,  fion0mer,u  ™pt,u  ki,i  i  =  1-4 
Subject  to  0  <C  £rj  <  1,  0  <  ^ionomer, i  <  1,  0  <C  £solid,i  ^  1  i  =  1— 4 

The  volume  fractions  of  the  voids,  ionomer,  and  solids  in  the 
catalyst  layer  are  a  function  of  the  optimization  variables.  For  all 
the  CLs  (i  =  1-4) 

Hlpt,  i  1  —  fpt,  i  ,  fpt,i  ,  f ionomer, i 

£r,i  =  1-7 - t—  - -  +  —  +  tz - - f -  (1) 

tCL,iJPt,i  Pc  PPt  l  A  Jionomer,i\Pionomer 


Table  5 

Performance  comparison  between  base  case  and  optimized  designs. 


Voltage  (V) 

Current  density  (mAcm-2) 

Base  case  design 

Low  Pt 

High  Pt 

0.9 

12.6 

13.86 

23.65 

0.7 

457.3 

462.68 

548.47 

0.5 

1223.71 

1419.1 

1481.76 

o  _  1  Wlpti  f ionomer, i 

£  ionomer, i  —  I  ~  r  A  _  f 

L CL, iP ionomer  JPt,i  |_  1  Jionomer,i 

£  solids, i  =  1  ~  &r,i  ~  &  ionomer,  i 

3.2.  Discussion  of  the  results 

The  comparison  between  the  base  case  design  and  optimized 
conditions  at  various  operating  voltages  for  both  the  cases  is  shown 
in  Table  5.  In  Table  5,  it  is  observed  that  the  performance  enhance¬ 
ment  is  quite  less  at  high  operating  voltages  (low  current  densities) 
when  the  maximum  platinum  loading  is  the  same  as  that  of  the 
base  case  design.  On  the  other  hand,  the  enhancement  in  the  per¬ 
formance  is  almost  double  ( 1 2.6  mA  cm-2  to  23.65  mA  cm-2 )  when 
the  upper  bound  on  the  platinum  loading  is  high  (1  mg  cm-2)  as 
seen  in  the  first  row  of  Table  5.  This  is  due  to  considerable  reduc¬ 
tion  of  the  activation  losses  by  using  high  platinum  content  at  low 
voltages.  With  an  increase  in  the  current  density,  two  trends  are 
observed.  First,  for  the  same  low  platinum  loading  (0.4  mg  cm-2), 
the  optimized  design  improves  the  performance  by  roughly  16% 
( 1 223.7 1-1419.1  mA  cm-2 ).  Second,  the  impact  of  higher  platinum 
loading  decreases  as  the  operating  voltage  becomes  lower  mainly 
due  to  the  increase  in  the  ohmic  and  diffusion  losses  (high  Pt  case). 
The  results  also  show  that  the  optimum  platinum  loading  hit  the 
upper  bound  (0.4  mg  cm-2)  in  the  low  Pt  loading  case.  In  the  high 
Pt  case,  though  the  optimum  loading  reached  the  upper  bound  at 
low  current  densities,  it  decreased  at  high  current  densities.  These 
results  indicate  that  the  cell  performance  can  be  enhanced  for  the 
same  platinum  loading  with  the  optimized  design.  However,  since 
the  best  performance  can  be  achieved  at  higher  platinum  loading 
and  the  objective  of  the  first  study  is  to  maximize  cell  performance, 
detailed  analysis  presented  below  is  performed  with  high  Pt  loading 
case  only. 

The  comparison  between  the  base  case  design  and  optimized 
conditions  is  shown  in  Fig.  2.  Optimized  cell  current  densities  are 
shown  as  dashed  lines  as  each  current  density  is  obtained  using  a 
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Fig.  2.  Performance  comparison  between  base  case  design  and  optimized  condi¬ 
tions. 
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Fig.  3.  Optimum  distribution  of  (a)  weight  fraction  of  platinum  on  carbon,  (b)  weight 
fraction  of  ionomer  and  (c)  platinum  loading. 


different  set  of  catalyst  layer  parameters.  Here,  the  base  case  design 
represents  a  single  CL  of  thickness  20  pan.  For  the  base  case  design 
performance  curve,  the  design  parameters  are  uniform  in  the  entire 
catalyst  layer.  The  optimum  values  of  the  CL  parameters  at  various 
voltages  are  shown  in  Fig.  3.  In  these  plots,  CL1  corresponds  to  the 
MPL/CL  interface  and  CL4  represents  the  CL/MEM  interface. 

Fig.  2  shows  that  there  is  significant  improvement  in  the  cell 
performance  in  all  limiting  regions.  The  improvement  in  the  current 
density  due  to  optimization  is  found  to  be  about  33%  at  0.4  V.  At  low 
current  densities,  the  platinum  loading  in  the  optimal  design  hit 
the  upper  bound  in  all  the  CLs  in  order  to  minimize  the  activation 
losses.  A  similar  trend  is  observed  earlier  in  the  studies  with  single 
CL  [10]. 

As  the  ohmic  and  mass  transport  losses  become  significant  at 
medium  and  high  current  densities,  an  increase  in  the  ionomer  and 
void  fractions  and  reduction  in  the  platinum  loading  are  observed 
in  all  the  CLs.  The  optimized  design  variables  which  support  this 
observation  can  be  seen  from  Fig.  3.  An  increase  in  the  optimum 
weight  fraction  of  ionomer  can  be  observed  in  Fig.  3b  with  the 
increase  in  current  density.  Overall  optimum  platinum  loading 
reduces  by  16%  when  the  operating  voltage  is  reduced  from  0.6  V 
to  0.4  V  (Fig.  3c).  The  enhancement  in  the  cell  current  density  for 
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Fig.  4.  Optimum  distribution  of  void  fractions  in  the  CLs. 


the  optimal  design  is  more  significant  at  high  current  densities. 
This  is  due  to  a  reduction  in  the  cathode  overpotential  due  to  opti¬ 
mum  distribution  of  voids  and  ionomer  in  all  CLs.  For  example  at 
a  voltage  of  0.5  V,  the  optimum  void  and  ionomer  fractions  in  CL1 
through  CL4  are  0.29,  0.24,  0.19,  0.13  and  0.27,  0.37,  0.45,  0.56, 
respectively.  The  optimum  void  and  volume  fractions  of  all  the  CLs 
at  various  operating  voltages  are  shown  in  Figs.  4  and  5,  respec¬ 
tively.  As  the  flux  of  oxygen  decreases  from  the  MPL/CL  interface 
towards  the  membrane  due  to  the  electrochemical  reaction,  the 
optimum  void  fractions  also  decrease.  Since  more  amount  of  oxy¬ 
gen  is  required  at  higher  current  densities,  a  gradual  rise  in  the 
optimum  void  fractions  is  observed  from  0.7  V  to  0.4  V  in  all  the 
CLs. 

The  trend  of  the  optimal  ionomer  fraction  is  opposite  to  that 
of  the  void  fraction  due  to  the  consumption  of  protons  from  the 
CL/MEM  interface  towards  the  MPL/CL  interface.  The  optimal  dis¬ 
tribution  of  the  ionomer  fraction  is  shown  in  Fig.  5.  The  thicknesses 
of  the  individual  CLs  are  found  to  hit  the  upper  bound  of  5  [xm  in 
both  high  and  low  Pt  loading  cases.  This  study  shows  that  the  overall 
CL  thickness  of  an  optimized  multiple  CL  PEMFC  can  be  higher  than 
the  typical  CL  thickness  of  10  [xm  or  less  used  in  the  conventional 
single  CL  PEMFCs. 

Although  the  optimum  design  values  of  the  CLs  are  different 
for  each  operating  voltage,  an  optimized  MEA  can  be  prepared  only 
with  one  set  of  these  parameters.  To  choose  such  a  combination,  the 
optimum  design  parameters  corresponding  to  0.5  V  (power  den¬ 
sity  is  maximum  at  this  voltage  in  the  base  case  design)  have  been 
selected.  The  steady  state  simulations  are  carried  out  with  this  opti¬ 
mum  design  to  find  the  power  densities  throughout  the  operating 
range.  The  comparison  between  the  power-voltage  curves  of  the 
base  case  design  PEMFC  and  the  PEMFC  with  optimally  designed 
multiple  CLs  is  shown  in  Fig.  6.  The  profile  of  oxygen  partial  pres¬ 
sure  in  the  base  case  design  and  optimized  designs  is  shown  in 
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Fig.  5.  Optimimum  distribution  of  ionomer  volume  fractions  in  all  CLs. 
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Fig.  6.  Performance  curves  for  the  base  case  design  and  for  the  optimal  design  at 
0.5  V. 

Fig.  7.  The  partial  pressure  of  oxygen  is  higher  in  most  part  of  the 
CLs  in  the  optimized  multiple  CL  design.  As  explained  above,  this  is 
due  to  the  optimum  distribution  of  voids  in  the  multiple  CL  design. 

3.3.  Performance  comparison  between  the  optimized  single  CL 
cathode  and  multiple  CL  cathode 

Fig.  8  shows  a  comparison  between  the  performance  of  the  opti¬ 
mized  PEMFC  with  four  CLs  and  the  optimized  PEMFC  with  single 
CL.  In  Fig.  8,  the  vertical  axis  represents  the  percentage  increase 
in  the  cell  performance  compared  to  the  base  case  design.  The  blue 
and  red  colors  represent  the  performance  of  the  optimally  designed 
PEMFCs  with  single  CL  and  multiple  CLs,  respectively.  About  3-6% 
increase  in  the  cell  performance  is  observed  with  the  PEMFC  with 
optimally  designed  multiple  CLs  in  comparison  to  the  PEMFC  with 
optimally  designed  single  CL  throughout  the  polarization  range  for 
the  same  platinum  loading.  Fig.  9  shows  the  difference  between 
the  power  density  achieved  using  optimized  multiple  CL  PEMFC 
and  optimized  single  CL  PEMFC  in  the  entire  polarization  range. 
As  seen  in  the  figure,  the  power  density  of  the  optimized  multiple 
CL  PEMFC  is  23  mW/cm2  higher  at  0.5  V.  For  a  large  stack  totaling 
thousands  of  cm2  of  area,  the  gain  in  total  power  generation  can  be 
very  high. 
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Fig.  8.  Percent  increase  of  current  density  of  optimzed  designs  of  single  CL  and 
multiple  CLs  over  base  case  design. 
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Fig.  9.  Power  density  enhancement  due  to  optimized  multiple  CL  design  over  opti¬ 
mized  single  CL  design. 


3.4.  Optimization  formulation  II 

In  this  formulation,  the  platinum  loading  is  minimized  while 
maintaining  the  same  performance  as  the  base  case  design.  In  addi¬ 
tion  to  the  inequality  constraints  in  formulation  I,  an  additional 
constraint  on  the  cell  performance  is  considered  ( iopt  =  hase )•  Here, 
iopt  is  the  cell  current  density  at  the  optimum  design  and  ibase  is  the 
base  case  design  current  density. 

Objective  function:  Maximization  of  iceiilmpt  at  a  given  VCeu 

Decision  variables  :  fptJ,  fonomerj,  ta,i  i  =  1-4 

Subject  toO  <  £r  i  <  1,  0  <  Sionomer,i  <  1,  0  <  SsoMA  <1  1  =  1-4,  lopt  =  ibase 

Optimization  is  performed  at  low  and  high  currents,  and  the 
optimum  platinum  loadings  are  shown  in  Fig.  10.  In  Fig.  10,  black 
and  blue  colors  represent  the  platinum  loading  of  the  base  case 
design  and  optimum  designs  with  four  CLs,  respectively.  There  is 
a  significant  increase  in  the  current  generation  per  mg  of  catalyst 
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Fig.  7.  Comparison  of  oxygen  partial  pressure  between  base  and  optimized  designs. 


Fig.  10.  Comparison  of  optimized  and  base  case  design  platinum  loadings. 
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Fig.  11.  Optimum  platinum  loading  with  single  and  multiple  CLs  cathode. 

loading.  In  other  words,  there  is  a  considerable  reduction  in  the 
platinum  loading  for  the  same  performance.  When  the  current  den¬ 
sity  is  0.15  A  cm-2,  the  optimum  platinum  loading  is  0.33  mg  cm-2 
against  the  base  case  design  platinum  loading  of  0.4  mg  cm-2.  This 
is  a  17%  reduction  in  the  platinum  loading.  When  the  current  den¬ 
sity  increases  further  to  0.73  A  cm-2,  the  platinum  reduction  is 
about  60%.  The  platinum  reduction  is  more  at  high  current  densities 
because  of  lower  activation  losses  and  stronger  effect  of  ohmic  and 
mass  transport  losses.  Once  again,  the  thicknesses  of  the  individual 
CLs  are  found  to  hit  the  upper  bound  of  5  p,m. 

The  obtained  results  are  compared  with  the  optimum  design 
of  a  single  CL.  The  red  color  in  Fig.  10  represents  the  optimum 
platinum  loading  required  to  produce  the  same  current  density 
(0.1 5  A  cm_2and  0.73  A  cm-2 )  using  the  optimized  single  CL  design. 

Fig.  1 1  shows  the  optimized  platinum  loadings  of  single  CL  and 
multiple  CL  PEMFCs  at  0.8  V  in  comparison  to  the  base  case  design. 
The  base  case  design  conditions  are  given  in  Table  3.  Here,  opti¬ 
mization  is  carried  out  with  the  steady  state  model  of  a  PEMFC 
cathode  containing  one,  two,  three,  and  four  catalyst  layers.  The 
total  thickness  of  the  reaction  medium  is  kept  constant  at  20  |jim. 
The  thickness  of  each  CL  is  calculated  such  that  the  overall  thickness 
is  20  |jim.  For  example,  while  considering  two  CLs,  the  thickness 
of  each  CL  is  10p,m.  When  a  single  CL  (20p,m)  is  considered,  a 
4%  reduction  in  platinum  loading  is  achieved.  With  two  CLs  (each 
of  10|jim  thickness),  the  reduction  is  9%.  The  reductions  in  plat¬ 
inum  loading  with  3  and  4  CLs  are  1 2.5%  and  1 5.5%,  respectively.  As 
the  number  of  catalyst  layers  increases,  the  design  variables  also 


increase.  As  a  result,  flexibility  in  selecting  the  optimum  design 
variables  also  increases. 

4.  Conclusions 

In  this  paper,  two  optimization  studies  of  a  PEMFC  with  an  inno¬ 
vative  multiple  CL  cathode  using  a  two-dimensional  two-phase 
model  are  presented.  The  decision  variables  used  in  these  studies 
are  the  design  parameters  of  all  CLs.  In  the  first  optimization  study, 
the  cell  performance  is  maximized.  The  optimization  improves  the 
cell  current  density  by  about  1 5%  in  the  high  current  density  region 
and  by  about  85%  in  the  low  current  density  region.  The  PEMFC  with 
optimized  multiple  CLs  shows  an  improvement  of  3-6%  in  the  cell 
performance  over  the  optimized  PEMFC  with  single  CL.  In  the  sec¬ 
ond  optimization  study,  the  platinum  loading  has  been  minimized 
for  a  given  performance.  The  optimization  results  in  a  reduction  of 
17-60%  platinum  loading  at  various  current  densities  in  compar¬ 
ison  to  the  base  case  design  without  sacrificing  the  performance. 
When  the  optimization  is  performed  with  single,  two,  three,  and 
four  catalyst  layers  at  0.8  V,  it  is  observed  that  a  higher  reduction 
in  the  platinum  loading  can  be  achieved  with  an  increase  in  the 
number  of  CLs.  This  study  shows  that  an  optimized  multiple  lay¬ 
ered  PEMFC  is  a  promising  option  to  improve  the  performance  of 
the  PEMFCs  and  to  decrease  its  cost. 
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